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ABSTRACT 



Context. Alfven waves and ion-cyclotron absorption of high-frequency waves are frequently brought into models devoted to coronal heating 
and fast solar-wind acceleration. Signatures of ion-cyclotron resonance have already been observed in situ in the solar wind (HELIOS 
spacecrafts) and, recently, in the upper corona (UVCS/SOHO remote-sensing results). 

Aims. We propose a method to constrain both the Alfven wave amplitude and the preferential heating induced by ion-cyclotron resonance, 
above a partially developed polar coronal hole observed with the SUMER/SOHO spectrometer. 

Methods. The instrumental stray light contribution is first substracted from the spectra. By supposing that the non-thermal velocity is related to 
the Alfven wave amplitude, it is constrained through a density diagnostic and the gradient of the width of the Mg x 625 A line. The temperatures 
of several coronal ions, as functions of the distance above the limb, are then determined by substracting the non-thermal component to the 
observed line widths. 

Results. The effect of stray light explains the apparent decrease with height in the width of several spectral lines, this decrease usually starting 
about 0. 1-0.2 Rq above the limb. This result rules out any direct evidence of damping of the Alfven waves, often suggested by other authors. 
We also find that the ions with the smallest charge-to-mass ratios are the hottest ones at a fixed altitude and that they are subject to a stronger 
heating, as compared to the others, between 57" and 102" above the limb. This constitutes a serious clue to ion-cyclotron preferential heating. 



Key words. Sun: corona - Sun: UV radiation - Line: profiles 



1. Introduction 

Ion-cyclotron resonance is often suggested as one of the pos- 
sible processes that accelerate the fast solar wind and heat 
the corona. Different modelling approac hes are used, like hy- 
brid simu lations (e.g. Liewer et al.lll999tlOfmanl2004l). kinetic 
models jl senberg et alj l2Q0lk llsenbergl |200lt iMarsch & Tul 
32), c 



the different local ion gyrofrequencies decrease with altitude 
and match the wave frequency, taking into account the Doppler 
effect due to the particle velocity ( for the theoretical aspects o f 



the cyclotron absorption, see e.g. iHollweg & Isenbe rg 2002). 
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Vocks & Marsch 2002), or fluid models. In this last cate- 
gory ICranmer et alj dl999al) for example, include several ion 
species and show that it is necessary to regenerate the cy- 
clotron waves through the wind to provide enough energy to the 
wind. This may be done, for example, by includ ing an Alfven 
w ave turbulent cas cade, as in lTu & Ma rsch ( 1997) (two fluids) 



Hu et alj j2OQ0h (fou r fluids ), while iMarkovskiil d200ll) and 



Markovs kii & Hollwegl (2002) study the direct generation of 
cyclotron waves via plasma instabilities. 

In the models that use turbulent cascade, the Alfven wave 
power is distributed into a power law spectrum, with ion- 
cyclotron waves being its high frequency limit. On the one 
hand, these waves are gradually absorbed by the particles as 
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On the other hand, they replenish through the energy cascade 
process from low to high frequencies. The large decrease in 
the magnetic field intensity because of the rapid expansion of 
the flux tubes can favour the resonance occurring at a relatively 
low altitud e above the solar surfac e, a s in th e so-called "fun- 
nels" (e.g. lHackenberg et al. 2000;O|200l). In this respect, 
the superradial expan sion of the flux tubes in coronal holes 
( Kop p & Holzet Il976l) . which are thought to be the source of 
the fast solar wind, is also important. 

In some models, Alfven waves are dire ctly di s sipate d 
to provide energy to the wind. According to IParken (1199 II) . 
though, this dissipation is not significant below 5 R , so that 
it is more likely that Alfven waves stay undamped in the low 
corona. 

Several signatures of wave-particle interaction, which can 
be attributed to the ion-cyclotron resonance, have already been 
observed in the solar wind. The in situ observations of the 
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Helios spacecrafts, for example, showed differential speed of 
minor heavy ions, anisotropies of temperatures (T , /T\\), and 
prefe rential heating of some ion species dSchwenn & Marsch 
199 lb . In the upper corona, the UVCS spectro-coronometer 
on board SOHO made it possible to observe the same kinds 
of signatures from 2.5 to more than 5 R Q from disc cen- 
tre, th rough spectro s copic studies coupled to some mod- 
elling dLi et al.lll998t Eohi et alJll998t Icranmer et al.lll999bl 
Esser et al.ll 19991) . In the lower corona, in contrast, there is lack 
of conclusive studies, and yet it is particularly important to 
verify whether any process that accelerates the fast solar wind 
starts acting in the lower corona. Observing this region is also 
important for constraining the energy actually dissipated by the 
cyclotron resonance and for differentiating the different theo- 
ries that have been proposed. 

In this article, we concentrate on the possibility of reveal- 
ing the preferential heating of the ion species having the lowest 
charge-to-mass ratios {q/m). Retrieving the temperature of dif- 
ferent ion species can be made via the analysis of the widths of 
emission lines. It is worthwhile noting that the temperature di- 
agnostics based on line ratios are sensitive to the electron tem- 
perature and cannot help in measuring the ion one. 

Our paper is organised as follows. In Sect. [2] we present 
the usual interpretation of the coronal line widths and of their 
variations above the limb, as well as solutions proposed in the 
literature for distinguishing the thermal and the non-thermal 
components. In Sect. [3] we present the observations and the as- 
sociated data processing, with emphasis on the correction for 
the stray light contribution in the spectra. In Sect. @] we ex- 
plain the method we used to separate thermal and non-thermal 
Doppler effects in the measured line widths. The results are 
presented in Sect. [5] while Sect.[6]gives the conclusions. 

2. Retrieving the information contained in coronal 
line widths 

2.1. Interpretation of the VUV line widths of coronal 
ions 

The plasma observed above the limb is optically thin for most 
of the transition region and coronal lin es in the VUV range 
observed by instruments like SU MER dDere & Masonl 1 1993t 
Chae et al.lll998HLee et al.ll2000h . The two main contributions 
to the line width then come from thermal and non-thermal 
Doppler effects. The non-thermal velocity, also known as "un- 
resolved" velocity, £ was originally introduced to interpret the 
width of the transition region and coronal lines observed on 
the disc. In fact, there is an excess broadening as compared to 
the thermal component due to the formation temperature Tf of 
these lines. This additional contribution was first attributed to 



thermal contributions are Gaussian (ITu et al.lll998), the half- 



veloc ity fluctuations induced by acoustic waves (IBoland et al.l 
19731 before MHD waves or fluid turbulence were considered 



(e.g.lBoland et alJll975tlHassler et al.lll99oT:lDovle et ai1ll999t 



Harri son et al.ll2002l) . The term "unresolved velocity" in fact 



refers to motions on scales smaller than the resolution scale. 
One should also take into account the motions integrated on the 
depth of the line of sight (LOS), and on scales smaller than the 
temporal resolution of the instrument. If both thermal and non- 



width cr at 1/ yfe follows the formula given by Dere & Masonl 
dl993h : 



A 2 2kT 

n 9 ^ 

2c z m 



■f) 



(i) 



where A is the wavelength, c the speed of light, k the Boltzmann 
constant. The non-thermal velocity f is proportional to the 
mean velocity fluctuations integrated on the observed solid an- 
gle and the exposure time, < 5v 2 > 1 ^ 2 . The coefficient of pro- 
portionality, relatively close to 1, depends on the degrees of 
freedom of the fluid motion (e.g. lTu et al.lll998 ). 

Variations in the non-thermal velocity have been inves- 
tigated by se v eral a ut hors: (i) centre-to -l imb variat i ons by 



Mariska etajl (1 19781) . IChaeetalJ (1 1998 ) iLeeetalJ J2000b . 



Doyle et al. ~d2000 ). see also [krdelyi et~ 



(1998) for the inter- 



pretation in terms of MHD modes; (ii) variation depending on 
the nature of the observed region (coronal holes, "quiet Sun", 
active region); and (iii) radial variation out of the limb, which 
is of interest for this article and which we will present in the 
next section. Depending on all these different conditions, the 
value of £ varies roughly from 10 to 30 km s _1 . 

From on-disc observations, we retain the results of 



Chae et al.l (11998b . which assumed that the temperature of each 
ion species that they observed was the corresponding forma- 
tion temperature (Tf). They consequently deduced the non- 
thermal velocity from the line widths, and plotted it as a func- 
tion of Tf . The same kin d o f analysis was also per formed by, 
e. gjBoland et alJd 1975b and lDere & Masonl dl993b . The curve 
of IChae et al.l dl998b increases from » 5 kms~' at 10 4 K up 
to as 30 kms _I around 3 x 10 5 K (Ovi), then decreases to 
ss 17 kms~' at 10 6 K (Fexn). One can roughly interpret this 
curve in terms of LOS integration: different Tf of ions observed 
at the same pointing on the disc can be associated to differ- 
ent altitudes above the solar surface. The curve then traces the 
value of the non-thermal velocity thr ough the different t emper- 
ature layers of the solar atmosphere ( Mariska et al.lll978 >. 



2.2. Past observations of line widths above the solar 
limb 

Most authors who have observed above the limb, in coronal 
holes, report that the line widths first increase with the altitude 
before presenting a plateau or, s ometimes, starting to decrease 
at 1.1 R Q from disc centre: e.g. lHassler et al. (1990), with the 



Mgx 609 and 625 A lines up to 1.2 R G , or iBaneriee et al 



dl998h . with SUMER/SOHO observing the Si vm 1445 A line. 
These authors interpreted this behaviour as caused by Alfven 
waves, with the increase due to flux c onservation and decr ease 
due to damping. Using CDS/SOHO, lO'Shea et al] d2003b ob- 
served the width of the Mgx 625 A line starting to decrease 
around 1.06 R Q , the exact position varying with the data 
sets, but alw a ys for the same value of the line widths. But 
O'Shea et al.l d2004l) . when observing the Si xn 520 A line in 
addition to the two Mgx 609 and 625 A lines, found a dis- 
crepancy between both species and concluded that the be- 
haviour of the Mg x lines results from photoexcitation grow- 
ing as the altitude increases. Still above a coronal hole, using 
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SUMER, [Lee et all d2000l) found a link between the altitude 
where the line widths change their behaviour and the intensity 
scale height. Nevertheless, they used a lot of particularly cold 
lines, which they could not observe more than 20 " above the 
limb. They also found a curve E, — /(7f) as did IChae et all 
dl998h . 



Above a quiet Sun region at the equator, lHarrison et al 



(120021) observed a decrease in the line width of the Mgx 
625 A line, with CDS, starting at about 1.04 Rq, the width 
being narrower than above coronal holes. They linked this be- 
haviour to the dampi ng of Alfven waves in close field loops. 
Wilhelm etall d2004l) . using SUMER still at the equator but at 
a different date, did not find the same behaviour. The results 
of both instruments were reconciled after joint observations 
( Wilhelm et al.ll2005[) that revealed no decrease, but a slight 
plateau. Besides. IPeter & Vocksl (12003) observed a local maxi- 
mum of the widths of transition region lines (O v and S vi), less 
than 10" (~ 0.01 R Q ) above the limb, which they attribute to the 
dissipation of ion-cyclotron waves, the decrease at higher alti- 
tudes being due to collisional cooling when the ion-cyclotron 
dissipation vanishes. 

Singh et ail d2003albl) analysed ground-based observations 



of Fe x, Fe xm, and Fe xiv lines (visible and infrared), up to 
140" in streamers. On average, the Fex 6374 A line width in- 
creases with altitude, those of the Fexm 10747 and 10798 A 
lines increase sligthly less, whereas that of the Fe xiv 5303 A 
line decreases slightly. Concluding that the line-width be- 
haviour depends on the formation temperature, they reject the 
interpretation in terms of variation of Alfven wave amplitude 
and prefer to explain it by an LOS effect in a thermally s trat- 
ified corona. Still from the ground, IContesse et al.1 ((2004) ob- 



serve no clear decrease in the width of the Fe xiv 5303 A line, 
but instead a plateau up to 1 .5 R Q at the equator (east and west) 
and up to 1.2 R Q above both poles, the observed line widths 
being larger above the poles than above the equatorial regions. 

Using UVCS/SOHO, it is possible to observe several so- 
lar radii above the solar limb, essentially the H 1 Ly a and the 
Ovi lines at 1032 and 1037 A. The observations showed that 
widths generally increase more with the altitude in the upper 
part of the cor ona, i.e. above 2 R from disc centre, than in 
the lower one dLi et alJll998l: iKohl et alJll998t ICranmer et alJ 
1999bHEsser et al|ll999l:lFrazin et alJll999l:lDovle et alJll999t 
Morgan & Hab bal2004l) . 



2.3. Distinguishing thermal and non-thermal Doppler 
effects 

The line width contains information about the result of poten- 
tial heatings (the temperature T x of the different ions species) 
and about the possible source of heating, that is to say, 
the Alfven wave amplitude via direct dissipation or via ion- 
cyclotron resonance after a turbulent cascade. All difficulties 
come from these two pieces of information merged in only one 
observable. Most methods used in the literature consist in mak- 
ing an assumption about one of the above quantities, depending 
on the purpose of the authors, i.e. constraining the Alfven wave 
amplitude or revealing any preferential heating. 



One can find mostly two kinds of hypotheses on the ion 
temperature: the isothermal or the formation temperature hy- 
poth eses. One example of i sothermal hypothesis can be found 
in iFeldman et alj J2000l) or boschek et all d200lb . They deter- 



mine the electron temperature through the emission measure 
and assume that all species are in thermal equilibrium. Then, 
they find different values of £ for the different ion species and 
discuss about possible errors to e xplain why th ey do not ob- 
tain a unique value for In fact, iMoranl d2003l) shows for a 
coronal hole that it is not possible for all ion species to have 
the same temperature and the same non-thermal velocity at the 
same time. 

Numerous authors prefer to assume that the temperature 
7^ of each ion is t h e formation temperat ure of the line (e.g. 
Hass ler et al. 1990t Contesse et al. 2004) . But this supposes 



ionization equilibrium, while several authors think that the 
plasma above coronal hol es is weakly collisional dTu et alJ 
1998; Banerjee et al. 19981) . When a preferential heating pro- 



cess, like ion-cyclotron resonance, is at work on time scales 
smaller than the collisional cooling (cf. Table [TJ, the assump- 
tio n Tj - Tf no l onger makes sense. 



Dovl e et al 



(1999), for a coronal hole, use a more elab- 
orated hypothesis for the temperature. They suppose that the 
electron temperature Tq is equal to the formation temperature 
of the observed Si vm line and deduce the proton temperature 
7p from hydrostatic equilibrium, constrained by a density di- 
agnostic. Then, they suppose that 7\;j vm = Tp (and Tp > Tq) 
before they deduce £ from the Si vm line width. 

The first hypothesis about the non-thermal velocity ^ that 
can be made is that all ion species have the same ^. This is rea- 
sonable if this velocity comes from fluid motions, i.e. regard- 
less of the species; t his may be false if some sp ecies have high 
differential speeds (Hol lweg & Isenbergll2002l) . which should 
not be the case in the lo wer corona. The problem is then to 
determine the value of ISeely et alJ d 19971) adjust £ so that 
the temperature of the Fe ion species, i.e. the heaviest element 
they observe, is similar to that of ion species of neighbouring 
Tf, as they suppose that all ion species that have a common 
Tf are in thermal equilibrium. They then deduce the temper- 
ature of all the other ion species, which appear to be larger 
than their own formation temperature; the lower the forma- 
tion temperature, the larger the difference. These authors con- 
clude that this situation may come from out-of-equilibrium pro- 
cesses (heating or cooling phase) or from preferential heating. 
They also conclude that the u nderestimation of the ion temper- 
atures led Hassler et al. (1990) to overestimate £ at altitudes of 
100-200arcsec. 

When loo king for the si gnature of preferential heating in 
coronal holes, Tu et alJ dl998l) cautiously give an upper and a 
lower limit to The lower one being 0, they set the upper one 
by setting to zero the thermal width of the lines from the heav- 
iest ion having the smallest Doppler width (Fe xn) in their data 
set. Then, they constrain the temperature of all ion species, and 
note a stronger heating of ions having the smallest q/m between 
about 1.04 and 1.18 R . 

One possibility for better constraining £ is to suppose that 
the non-thermal velocity is solely due to undamped Alfven 
waves. To conserve the wave energy in an expanding flux tube, 
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Table 1. Comparison of the characteristic cyclotron time of 
some species with their heat exchange time with protons, in 
the lower corona. 



Species 


q/m* 


To (s) f 




T K (S) c 


e~ 


1758 


4x 10~ 


-7 


120 


P + 


1 


6x 10 


-4 


3 


Mgx 


0.37 


2x 10- 


-3 


1 


Fex 


0.16 


4x 10- 


-3 


2 



Charge-to-mass ratio, normalized to the proton one. Mg x and Fe x 
are the ion species having the highest and lowest values of q/m, 
respectively, in the spectrum observed by SUMER. 
t c = 2n-^, with q the charge and m the mass of the particle, and 

B the magnetic field. 

Heat exchange time with protons Tg p and t? (e.g. IWessorj[2004l . 
p. 68). 

All these times are calculated with the following conditions: B = 
1 gauss, electron density n e = 10 14 rrr 3 , and proton and electron 
temperatures both equal to 10 6 K (/3 ~ 0.35). 




while the mass dens ity p decreas es, the wave amplitude must 
increase so that (e.g. lMoranll200ll) 



-1/4 



(2) 



The line-width variation wi th the altitude has been analysed 
with respect to this ef fect inlHassler et alJ I l990) for Mg x and 
Baneriee et ail dl998l) and lPekiinlu et all d2002l) for Si vm. 

Conclusions that can be drawn from the study of VUV so- 
lar line widths appear to strongly depend on the method used 
to separate the thermal from the non-thermal contribution. Two 
points emerge from past studies: (i) the probable existence of 
a non-thermal velocity whose radial behaviour above the limb 
mimics that of undamped Alfven wave amplitude. The ques- 
tion of whether these waves begin to damp at a certain point 
above the solar surface does not seem to have been solved, (ii) 
Another process may be at work that affects different lines in 
different ways. A q/m dependence is suspected that could be 
associated to ion-cyclotron resonance. 

The purpose of this article is to verify and constrain both 
of these two points without using hypotheses that a priori pre- 
clude one of them. We will focus on observations of a polar 
coronal hole, below 1.2 R Q . To avoid bias resulting from the 
line formation (e.g. photoexcitation) and from unknown line 
blendings, it seems useful to observe as many lines as possible. 
A wide range of masses and q/m of the ion species also favour 
the distinction between the relative contribution of thermal and 
non-thermal velocity in the observed line widths. 

3. Observations and data processing 

The VUV Spectromete r SUMER on board th e SOHO space- 
craft is described in Wilhelm et al. (1 1995b . IWilhelm et al 



(119971) . and lLemaireetalJ (U997I) . Table |2] briefly presents all 
data sets used in this article. The main one, called set 1, consists 
of ten spectral ranges, observed using detector A on 30 May 
2002 during MEDOC Campaign #9. The lines of interest are 
listed in Table[3] The complete observation was about 13 hours 



Fig.l. EIT/SOHO context image for data set 1 (195 A chan- 
nel). The spatial coverage of the SUMER slit is shown. 

Table 2. Description of data sets. 



Data set a 


Date 


Position b 


R G (arcsec) c 


1 


30 May 2002 


X=0, Y=1150 


956.5 


2 


28 May 2002 


X=0, Y=-1150 


956.5 


3 


19 Nov 2003 


X=0, Y=1145 


980.0 


4 


12 June 2004 


X=0, Y=-1140 


960.0 



a Observations are all made with the 1" x 300" slit, and detector A. 

In data sets 2, 3, and 4, only the Mg x 625 A line was analysed in 

the present paper, for comparison purpose. 
b Position of the centre of the slit, in arcseconds from disc centre. 
c Solar radii as seen from SOHO. 



long. The centre of the slit was pointed at X=0", Y=l 150" (i.e. 
at the North Pole). The 1" x 300" slit was used, oriented in the 
north-south direction, covering a wide range of altitudes at the 
same time (from about 50" to 350" above the solar limb, with 
an apparent solar radius of 956"). The width of 1" reaches a 
compromise between the photon statistics and the instrumental 
width. 

The EIT/SOHO images on that day (Fig. |T|l show that the 
coronal hole appearing on the solar disc at the North Pole is 
not developed well. Inspection of EIT movies of the preceding 
and following solar rotations reveals that the coronal hole is far 
from being symmetric, and highly elongated along one partic- 
ular direction towards the equator. Nonetheless, it appears that 
some part of the coronal hole is lying right under the position 
we are pointing at in the plane of the sky, and LASCO images 
show no streamer on the north-south axis. However, we cannot 
dismiss the possibility of a contamination of the LOS by any 
plasma belonging to non-hole (i.e. "Quiet Sun", QS) part of the 
corona. To gain more confidence, we compared the behaviour 
of the width of the Mg x 625 A line (observed at 1249.88 A in 
second order) in three other data sets. Data set 2 corresponds to 
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the south pole observed two days before data set 1 (same alti- 
tude range); the EIT survey shows that there is no coronal hole 
at the south pole at this period. Data sets 3 and 4 correspond 
to other polar coronal holes (CH). Figure [2] shows that the be- 
haviour of the line width in set 1 compares well to those in sets 
3 and 4, because it keeps increasing at low altitudes and then 
decreases at the highest altitudes. In contrast, data set 2 dis- 
plays a different behaviour, with a large plateau starting from 
60", and it wi ll be analysed more th oroughly in a future paper. 
According to Iwilhelm et all d2004 . the Mg x 609 and 625 A 
line widths are larger in coronal holes than in the Quiet corona 
(by more than 30% at 100"). Data sets 1 to 4 do not show such 
large differences. The maximum one is between sets 1 and 2, 
and the Mg x line width in set 1 is only larger than that of set 2 
by * 4% at 100". In fact, from 60" to 1 00", our values are 
higher than those of IWilhelm et al.1 ( 120041) in the Quiet corona 
by about 15% and lower than those in coronal holes by about 
15%. Nevertheless, t hey are closer to what was observed by 
Dosch ek et al. too lh in coronal holes, although a contamina- 
tion by non-hole material cannot be dismissed, especially be- 
cause the magnetic topology is far from being as simple as dur- 
ing solar minima. The density, discussed in Sect. 14. 11 is similar 
to what has been previously measured in coronal holes. 

A density diagnostic has been made using additional data 
acquired on 30 May 2002. The lines of Sivm at 1440 and 
1445 A were observed by SUMER, covering the same field 
of view as data set 1, with a 4" x 300" slit. This wider slit thus 
provides better photon statistics for the line ratio, at the expense 
of an optimal instrumental width. The total exposure time was 
about 2 hours and 50 minutes. 

Raw data were corrected using the standard SUMER proce- 
dures provided by Solar Software (flatfield and geometric dis- 
torsion of the detector). Note that the observation of faint lines 
with long exposure times is only possible b ecause the detecto r 
noise is low, « 3 x 10~ 5 countss" 1 px" 1 ( Wilhelm et alj[l997 ). 
which makes a; 0.1 counts px -1 for a 1 -hour exposure time. The 
final spectra were obtained through the addition of consecutive 
300 second exposures, after visual inspection of each exposure. 
The complete detector image was recovered (1024 x 360 pix- 
els, wavelength x spatial dimension Y). The instrument is stig- 
matic in the spatial dimension Y, with a theoretical spatial res- 
olution around 1". For better statistics, several spatial pixels 
were added before analysing the spectral profiles. These inter- 
vals appeared to be narrow enough to avoid any significant arti- 
ficial broadening of the lines owing to the residual distortion of 
the detector. A rough estimate of the altitude corresponding to 
each spectrum was provided through the pointing coordinates 
of the instrument. 



3.1. Correction from the instrumental stray light 

As SUMER is not a coronagraph (no occultor is included), the 
solar disc produces stray light during off-limb observations, es- 
sentially via the point-spread function properties induced by 
defects in the reflective surfaces, of spatial sizes of the same or- 
der of magnitude as the observed wavelength. Note that, what- 
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Fig. 2. Comparison of the line width (top panel) and the to- 
tal counts per second (bottom panel) measured for the Mgx 
625 A line (observed in second order; the measured width is 
then twice as large as the width that would be measured in 
first order), for the four data sets acquired above the solar poles 
(Table |2}. The error bars are only plotted for both extremities 
in the abscissas, for better visibility. The stray light contribu- 
tion has not been removed yet. The total counts in the line are 
discussed in Sect. 15. II 

ever the pointing is, the complete image of the Sun is perma- 
nently formed by the main parabola. 

This instrumental stray light is negligible for on-disc ob- 
servations, but the contamination increases as the altitude of 
observation above the limb increases, until it is responsible 
for the quasi-comple te photon rate detected above 1.6 R 
( Feld man et al. Il999h . This problem is even more critical in 
coronal holes than above the "quiet" corona, as the emissivity 
is lower. Consequently, each off-limb spectrum is the sum of 
the real spectrum emitted by the corona and of the stray-light 
spectrum. At each point of the image of the Sun produced by 
the mirror, the intensity of the stray light results from contri- 
butions convolved over the entire source (disc plus corona seen 
above the limb). But the main contribution comes from the disc, 
where the different regions (active regions, coronal holes, and 
"Quiet Sun") emit different spectra. This results in a stray-light 
spectrum that varies with the position and date of observations. 
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Fig. 3. Examples of stray-light spectrum predictions for a spec- 
tral range covering a; 40 A around the Si vm line at 944 A (data 
set 2). Each observed off-limb spectrum (in grey) is the sum of 
an instrumental stray-light spectrum (predicted in black) and of 
the real coronal spectrum. As the altitude increases (here, from 
60 to 134 and 253 arcsec above the limb), several cold lines 
become more prominent in the off-limb spectrum (like H i and 
S vi lines, here). Meanwhile, some coronal lines become more 
and more contaminated by the blending with neighbouring cold 
lines (like Si vm blended with S vi); blending can also be caused 
by the same coronal line when this one is prominent in on-disc 
spectra (e.g., for Mg x around 1250 A in second order). In these 
figures, spectra are in counts per wavelength pixels (with inte- 
gration over wider spatial ranges for increasing altitude). 



Every line that can be observed on the disc contributes to 
the stray-light spectrum observed above the limb, so that coro- 
nal lines (Mg x 625 A, Fe xn 1242 A, etc.) are not only blended 
with their own stray light, but they can also be blended with 
colder lines of neighbouring wavelengths. This last effect is 
usually underestimated. 

To correct for the instru mental stray light , we used a 



metho d alre ady de s cribed in IDolla et alJ (I2003h . IDolla et al 



d2004l) . and iDollal d2006l) . which enabled us to predict the 
stray-light spectrum contaminating the coronal spectrum ob 
served at any altitude. This predicted spectru m is then sub 



stracted from the off-lim b spectrum (see also Ban erjee et al 
1200(1 ISingh etaDl2003allbh . For the observed spectral ranges, 
a reference spectrum was taken well above the limb, so that 
the spectrum can be considered as purely due to stray light. All 
reference spectra were acquired with a long exposure (several 
hours), usually pointing at (X=0", Y=1700")- The same slit 
was used as for the off-limb spectra. The stray-light spectrum 
is then predicted at any lower altitude by determining the fac- 
tor that makes the total counts in a reference cold line coincide 
in both the observed off-limb spectrum and the predicted stray- 
light spectrum. If no line offers sufficient statistics, one can use 
the continuum, such as for the Fe x/Fe xi spectral range, near 
1465 A. Some examples of stray light prediction ar e shown in 
Fig.[3](see also lDolla et d]|2003il2004lDollall2006l) . 

We showed in IDolla et alJ (12003b that reference spectra 
taken directly on the disc pointing at a "quiet Sun" region, and 
averaged over all the spatial pixels observed by the slit, do not 
allow a precise prediction of the stray-light spectrum, because 
these spectra do not take every kind of region on the disc into 
account. With high-altitude spectra, on the contrary, the con- 
volution on the entire disc h as already been perf ormed. This 
also avoids modeling it like in lGabriel et alJ (120031) . Of course, 
the extrapolation of the high-altitude spectrum as a stray-light 
reference at any altitude is not totally correct, as the different 
contributions from the different areas of the disc are modified. 
This affects the ratios between the predicted lines, but the rela- 
tive error on t his ratio becomes more negligible as the altitud e 
increases (cf. lFeldman et al 1 19991: iDoschek & FeldmanE oOO). 
precisely as the absolute correction increases. 

During the observations, a high-altitude reference spectrum 
for the stray light was not taken for every spectral range ob- 
served at lower altitudes in the corona. In particular, judging 
from previous SUMER data and from spectral lines atlases, 
we thought that this was not necessary for the Si vm 1440 and 
1445 A lines. Nevertheless, Fig. [4] shows a bump in the blue 
wing, growing as the altitude increases. The same bump can 
be found in different data sets, so that statistical error is ex- 
cluded. A possible st ray light contamination, also noticed by 



Banerjee et al. (1998), may explain it. We suspect that another 



bump is present on the red wing. In fact, these bumps do not 
appear when using the 4" slit, probably because the wider in- 
strumental profil e smoothes the emission pr ofiles. 

According to lEkberg & Feldmanld2003h . the Sixi 604.15 A 
line is blended with an Fevn line. The blen ding line is a 



cold o ne (class "b", following the designation of iFeldman et al 



19971) . so that any direct blending is probably negligible above 



the limb. A stray light contamination is likely, but has been 
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Fig. 4. Si vm line profiles (1445 A) at 58" and 178" above a coronal hole (data set 1). No correction of the stray light has been 
applied, because no high-altitude reference spectrum was taken for that. The "bumps" marked by the arrows, getting larger with 
the altitude and present on other data sets, suggest that this line is slightly contaminated by one cold line or more, appearing as 
instrumental stray light. 



underestimated at the time of the data acquisition, and no high- 
altitude spectrum is available for the stray light determination. 



3.2. Gaussian width 

The half- width at 1/ -\[e is provided through a least-square 
Gaussian fit including a second-order polynomial. The spec- 
trum corrected for the stray light is in counts px~\ so as to use 
statistical errors as weights. The width in pixels is turned into 
angstroms by using the known dispersion of the spectrometer. 
The equivalent size of a pixel ranges from 42 to 45 mA in first 
order, depending on the observed wavelength. No precise cal- 
ibration of the li ne wavelength is nece ssary, and we used the 
values given by iFeldman et al. The fitting procedure 

provides a 1-sigma error on the width, which we set to 0. 1 pixel 
as a minimum. When the counts at line maximum were below 
150 or so, the fitting procedure became unstable, so that we 
only show results of lines presenting more than 150 counts at 
their maximum. 



4. Method of separating the two contributions in 
the coronal line widths 



The spectrum observed on the detector corresponds to the 
convolution of the spectrum arriving at the aperture of the 
spectrometer with the instrumental profile. This profile mainly 
comes from the finite width of the slit. Both off-limb and stray- 
light reference spectra are acquired with the same slit width, 
therefore the corrected s pectrum is s till convolved with the 
same instrumental profile (Dolla 2 0061) . 



The instrumental width, for detector A, does not depend 
on the wavelength. The procedure con_width_func t_3 in Solar 
Softw arfl updated following a study made by Chae et al. 
dl998l) . remove it through a deconvolution process, which 



proves to be equivalent to removing (quadratically) a Gaussian 
contribution of 42 mA (99 mA in FWHM). 



A newer procedure con_width_funct_4 differs only for detector 



As shown in lDolla et al J (120041) . that the non-thermal velocity is 
not well-constrained prevents direct conclusion about the pref- 
erential heating, by plotting either the temperature at a given 
altitude or the difference in temperature between two altitudes, 
as a function of q/m. In effect, owing to the small number of 
available lines, the ion species that are more likely to take ad- 
vantage of cyclotron heating, i.e. having low q/m ratios, are 
also those with low 1/ra ratios, that is to say, those whose width 
is more sensitive to the value of £ (cf. Eq.[T|i. 

Instead of making an assumption on the value of we pre- 
ferred to make an assumption on its nature. If we suppose that 
Alfven waves are responsible for the non-thermal velocity, then 
its radial variation can be constrained using Eq. [2] and a den- 
sity diagnostic. We also consider that the Mg x species, hav- 
ing the highest q/m ratio in our data set, is less likely to be 
heated through the cyclotron resonance, and we suppose that 
its line width will only increase because of an increase in the 
non-thermal velocity. This hypothesis will then be verified a 
posteriori on a given interval. 

4.1. Density diagnostic 

The Si vm spectra are produced by a running sum over 40 pix- 
els in the spatial direction of the detector. The total counts in 
the spectral ranges corresponding to both Si vm lines (1440 and 
1445 A) are integrated. The continuum contribution is then sub- 
stracted. The mass density p can be obtained through the elec- 
tron density ne, which is derived in Fig.|5]from the theoretical 
curve linking the line ratio (1445/1440) to the electron density 
dDoschek et alJI 19971 their Fig. 1). 

Above 270", the measured line ratio becomes lower than 
1, so that no solution is given by the diagnostic. The Sivm 
line ratio becomes far less sensitive at electron densities below 
10 7 cirT 3 and the 1445 A line suffer from the stray light men- 
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Fig. 5. Electron density ne with increasing altitude above the 
North polar coronal hole (data set 1). The dotted lines denote 
the error bars. 

tioned above. For all these reasons, we have less confidence in 
the results above 180". 



These density values are similar to those of iDoschek et al 



( 1997) measured in a polar coronal hole and at the same alti- 
tude range. They nevertheless decrease more rapidly with the 
altitude. 

4.2. Determination of £, at the lowest altitude 

Equation[T]can be rewritten as the quadratic sum of thermal and 
non-thermal velocity : 



2kT , 
m 



(3) 



By using the assumption that the Mg x 625 A line-width vari- 
ation is solely due to a variation in £ between positions ro and 
r (i.e. AT(r) = T(r) - T(ro) = 0, with ro the lowest position in 
data set 1 i.e. at the bottom of the slit), we can write 



A(v 2 (r)) 



v 2 (r)- 



The flux conservation of the Alfven waves yields 

m = (p(r)Y m 

£(ro) \p(r )j 

From Eq.[4]anf[5] we get 



A(v 2 (r)) 



( P(r) V 
Wo)/ 



-1/2 



1 



1/2 



(4) 



(5) 



(6) 



Thus, £(ro) can be determined for any r where the hypothesis 
Ar(r) = is true and where p(r) and v(r) can be measured. 

Figure [6] shows that the value of £(ro) is relatively constant 
for r in the range [70; 150] arcsec, which tends to prove that 
the AT = hypothesis is valid on that interval. As the error on 
£(ro) increases when the values of v 2 (r) and p(r) are close to 
that at ro, the lowest altitudes are associated with particularly 
large error bars. To get a higher precision, we made a weighted 



Fig. 6. Value of the non-thermal velocity £(ro) calculated at the 
lowest altitude observed in data set 1 (i.e. ho = 58" above the 
limb) for increasing altitude h — r - R Q (see Eq. |6]). We used 
the difference in width between r and ro, for the Mg x 625 A 
line. The dotted lines denote the error bars. 




100 150 200 

Altitude obove the limb (arcsec) 



250 



Fig. 7. Non-thermal velocity £ with increasing altitude above 
the North polar coronal hole (data set 1), derived by using the 
flux conservation of Alfven waves (Eq. [5J. The dotted lines de- 
note the error bars. 



average of £(ro) over the [90; 150] interval, reduced not to un- 
necessarily broaden the error bars. We got £(ro) = 15+2 kms -1 . 



4.3. Determination of £ and T at every altitude 

Once we have £(ro), £(r) comes from Eq. [5] (Fig. [7]), every 
where we know p(r). The Mgx 625 A line width and the 
AT = hypothesis were only used to determine the constant 
^(r ( >). The determination of £(r) now only requires the density 
diagnostic. Nevertheless, Eq.[5]is only valid if we assume that 
Alfven waves undergo no damping between ro and r and, of 
course, if Alfven waves are entirely responsible for the appear- 
ance of the non-thermal velocity. 
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Above 150", the slope of the density curve changes. This 
may be due to the presence of a different structure on the plane 
of the sky, e.g. a difference between plume and interplume. If 
so, one would expect that the flux conservation is no longer 
valid above 150", as the necessary continuity of the flux tubes 
is compromised. We nevertheless derived £ up to 270", keeping 
in mind that the results are less reliable above 150". 

Temperature for each ion species can now be determined 
from the observed linewiths (shown in Fig. [8]) by using Eq.Q] 
The altitude range is nonetheless restricted to positions where 
the non-thermal velocity is determined, i.e. below about 270". 
To produce each spectrum, we performed a running sum over 
40 spatial pixels. The results are shown in Fig.|9]and discussed 
below. 



5. Results and discussion 

5.1. No evidence of Alfven wave damping 

Before analysing more deeply the coronal line widths in terms 
of temperature and non-thermal velocity, we will discuss the ra- 
dial variation of the width as a whole. This is an important issue 
because, in past studies, the fact that the widths of coronal lines 
were decreasing above some given altitude was sometimes in- 
terpreted as the signature of the damping of Alfven waves. In 
fact, it appears that this decrease is essentially an effect of the 
instrumental stray light. 

Figure [8] shows some examples of the width of different 
lines, with increasing altitude above the limb. The Mg x 625 A 
and Si vm 1445 A lines, in particular, have been widely anal- 
ysed in the literature. When no stray light correction is made, 
the Mg x 625 A line presents the same behaviour as observed 
by previous authors: first, the line width increases, then reaches 
a maximum around 200" (~ 1 .2 R from disc centre), and starts 
to decrease. We find the same kind of behaviour for the Fe x 
1028 A or Fexn 1242 A lines. There are not enough statistics 
for the Si vm 1445 A line to extend the curve above 140". When 
we correct spectra from the stray light contribution in the case 
of the Mg x line, it is no longer possible to see any decrease, 
but maybe a flattening trend. For the Fe x and Fe xn lines, the 
stray light becomes so strong that there are not enough statis- 
tics to follow them as high as in the case with no correction; 
nevertheless, it appears an increasing trend for both lines. 

That the stray light contamination is responsible for reduc- 
ing the coronal line widths observed above the limb is well 
known, but appears to be underestimated in most studies with 
SUMER, as we find that the effect is already noticeable well 
below 200". One should question in particular any study that 
does not directly predicts the stray-light spectra. 

The explanation of this effect is rather simple: lines are al- 
ways narrower in stray-light spectra than in off-limb spectra; 
therefore, when the stray light contamination grows as the alti- 
tude increases, the observed line width is reduced. It is worth- 
while recalling that a stray-light spectrum is similar to any typ- 
ical spectrum that could be observed on the disc in the same 
wavelength range. Therefore two cases are possible: either the 
stray-light blending is caused by the same coronal or transition 
region line, when it is particularly prominent on the disc (e.g. 



for Mgx 625 A or Fexn 1242 A and 1349 A), or it is due to 
another line of neighbouring wavelength (e.g. Fe x lines that 
are blended with O i near 1028 A or with C i near 1463.5 A). In 
this case, the blending line is usually a colder and narrower one, 
dominant on the disc, while the hotter one is dominant above 
the limb. 

When no stray light correction is made, the observation of 
either a plateau or a decrease after the initial increase with the 
altitude, and the altitude where this is observed, is then sim- 
ply determined by the way the ratio of stray light to coronal 
emission evolves with the altitude. In general, lines having a 
lower formation temperature than the others present a radiance 
that decreases more rapidly with the altitude, and are conse- 
quently subject to more contamination. Thus, in Fig. [8] the 
uncorrected Fe x line width starts to decrease at a lower alti- 
tude (« 150") than that of Mgx (w 190"), which also starts 
to decrease at a lower altitude than that of Fexn 230"; cf. 
Fig.[l0]for the formation temperatures of these ions). But from 
one data set to another, the same line can also present differ- 
ent behaviours. One can notice in Fig. |2](upper panel) that the 
Mgx 625 A line width presents a plateau from 60" to 180", 
while it reaches a maximum at three different altitudes in data 
sets 1, 3, and 4. For these sets, the difference in altitude of the 
maximum present some correlation with the observed radiance 
(bottom panel). Note that if the observed radiance consists of 
the addition of the intrinsic and of the stray light radiance, the 
last one is negligible at low altitude, so that the curves give a 
good idea of the intrinsic radiance there. Thus, the curve starts 
to decrease before all the others in set 3 (~ 150") and reaches 
a plateau in set 4 (still around 150") before the maximum seen 
in set 1 (» 230"). This fact can be explained by an almost iden- 
tical stray light contribution in all data sets, dominating more 
rapidly in data sets where the intrinsic emission is smaller. 

There is, however, one reservation concerning the interpre- 
tation of the decrease in the line width in terms of stray light 
contamination: even after correction, there is still a decrease, 
or an almost constant value, for some lines in data set 1 : both 
Ovi, the Si xi 580.91 A, and the Sx 1196.20 A lines (Fig. [8] 
right panel). For these lines, there is no increase phase, as this 
decreasing trend begins at 57". Once the non-thermal velocity 
is removed, the derived temperatures decrease with increasing 
altitude (Fig. [9] third panel). We have no explanation for this 
behaviour, as these lines are not intense enough on-disc to be 
subject to photoexcitation, except for O vi. It is possible that 
these lines are affected by some blends. 

5.2. Signatures of preferential heating 

The ion temperatures in Fig. [9] are of the same order of mag- 
nitude as what is found in previous studies (e.g. lTu et al . 1998; 
Patsourakos et al.l2002h ; they all exceed 2 x 10 6 K at 57". They 
exhibit a wide range of values, which makes the isothermal 
hypothesis unlikely. The temperature of the three iron lines, 
for example, are arranged in the opposite order to their corre- 
sponding ionization temperature (Fig. |T0b . Then the LOS ef- 
fect associated to the stratification in temperature of the corona 
also appears unlikely. The temperature increases with the al- 
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Fig. 8. Line-width variation above the limb for some ions in data set 1. Spectra are created by a running sum over 40 spatial 
pixels. In grey, the value when the stray light is not corrected; the value after correction (in black) is not available for every ion 
(e.g. Si vm). The appearance of a decrease when some particular altitude is reached is explained by the stray light contamination 
(first panel, e.g. Fex blended with Oi, or Mgx, blended with the same line as stray light). The second panel shows some lines 
that start to decrease (or stay constant) from the lowest altitude that we observed (with or without correction), which apparently 
contradicts the behaviour of the other lines. For better visibility, the errors bars are only shown on both extremities of the curves, 
and the width of Fe x is shifted by -0.03 A (first panel), while that of S x is shifted by -0.02 A (second panel). 



titude for most of the lines, except in the bottom frame of 
Fig. [9] corresponding to the exceptions mentioned above. This 
implies some heating mechanism. The relative discrepancy be- 
tween both temperatures yielded by the Fe xn lines at 1242 and 
1349 A, below 130", may be explained by the blending of the 



1242 A line with an Si x line (cf. Table 1 in iFeldman et alj [Singhetal 



19971) . as well as by the error bars. The temperature of Fe x 
is very high, but it is confirmed at low altitude by a second Fe x 
line, at 1463 A (see Fig .P. 

To reveal a preferential heating, we plot in Fig. QT|the dif- 
ference in temperature as a function of the q/m ratio between 
57" and 102", for as many lines as possible in data set 1. The 
line widths and temperatures measured at these altitudes are 
presented in Table [3] It is important to note that, for lines seen 
in second order, the widths are twice as large as what would be 
measured in first order. This is taken into account when using 
Eq. [T]and the derived ones. We average over 30 pixels at 57" 
and 60 pixels at 102", to exploit some very faint lines. There 
were not enough statistics at 102" to analyse the Fex line at 
1463 A, so that we only keep the one at 1028 A. In the left 
panel, the non-thermal velocity has been taken into account, 
using results in Fig. [7] 15 and 17 kms~' at 57" and 102", re- 
spectively. The right panel corresponds to the case where no 
variation in non-thermal velocity is considered between 57" 
and 102", which is discussed in Sect. 15.31 The large error bars 
are due to the errors propagated from both the observed and 
the stray-light spectra to the line width (Gaussian-fit procedure) 
and the error on the determination of the non-thermal velocity. 
The weighted power law and linear fits show a trend for the 
ion species having the lowest q/m ratios to experience more 
heating than the others, which is compatible with ion-cyclotron 
preferential heating. Both kinds of fit curves are not very dif- 
ferent from each other. We do not intend to fit the data with any 
particular theoretical model. However, the power law tends to 
for high q/m, while the linear fit, with a slope of -1 .9 + 3.1 (in 



10 6 K per unit of normalized q/m), plummets towards -10 6 K 
for q/m = 1 (protons). Therefore, a power law seems to better 
represent a preferential heating of low q/m species, although 
the correct function is probably much more complicated. One 
can notice that the different line width behaviours observed by 

-i i ir-i 

(2003a b), which they link to the formation tem- 
perature (LOS effect), are also compatible with ion-cyclotron 
heating. It is worthwhile noting that they correct for the instru- 
mental stray light with the same method as we do. The case of 
Fe xiv, in their study, the line width of which decreases with the 
altitude, is nevertheless disturbing. But this is somehow similar 
to the behaviour of some lines in our data set 1 (Fig. [9] bot- 
tom panel). If it is a real variation in temperature, this apparent 
"cooling" process should be explained. We notice, though, that 
all these lines are associated to the highest q/m, meaning that 
they are less likely to be concerned by ion-cyclotron heating. 
Therefore, it is possible that all other ions experience the same 
kind of cooling, or a decrease in non-thermal velocity, but that 
this effect is counterbalanced by the preferential heating for the 
ions having the lowest q/m. 

A second kind of analysis consists in plotting the ion tem- 
peratures as functions of q/m, at a given altitude (Fig. 112b. It 
shows that the ion species having the lowest q/m have a higher 
temperature than the others. The trend is emphasised more at 
102" than at 57", which is consistent with the previous analy- 
sis. The particularly high temperature of Fe x tends to compress 
the graph, but the trend does not rely on this ion alone. To show 
that, we use a weighted linear fit including the Fe x data or not. 
Even if a straight line is obviously not the proper function to fit 
the data, it makes the comparison simpler, while no theoretical 
function is available. At 57", the slope is equal to -2.2 ±1.2 
with Fex and -0.9 +1.3 without it, still in 10 6 K per unit of 
normalized q/m. At 102", it gives -4.7 ± 1.7 with Fex and 
-2.7 ± 1.8 without it. 
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Fig. 9. Temperature with increasing altitude above the limb 
(data set 1 ; spectra are produced by a running sum over 40 spa- 
tial pixels). 



To interpret these curves, some points must be emphasised. 
First, the protons can be considered as a cooling source for all 
minor ion species that could be subject to any preferential heat- 
ing. They compose the largest part of the coronal plasma so that 
collisions of one given minor ion are much more frequent with 
protons than with other species. As the q/m ratio of the protons, 
equal to 1 with our normalization, is higher than those of all the 
minor ions, they are also least likely to experience cyclotron 
heating. Consequently, their temperature must be lower than or 



Fig. 10. Measured ion temperatures (data set 1, 57" above the 
limb) as fu nctions of the maximum ionization temperatures 
(taken from lArnaud & Rothenfluglll985 ). They are well above 
the corresponding ionization temperature (dotted line). 



equal to that of the ions having the highest q/m ratios, and the 
temperature of all other species is the result of a competition 
between preferential heating and cooling with the protons. By 
measuring differences in temperature, one only has access to 
the net heating. Second, most of the ion-cyclotron numerical 
simulations show that the species having the lowest q/m ab- 
sorb most of the wave power available, to the detriment of the 
others that cannot res onate efficiently yet (e.g. lCranmei 2000; 
Vocks & Marschl 2002 ). This could explain the high tempera- 
ture of Fe x as compared to the other species. 

That low q/m ions are already hotter than the others at 57" 
suggests that ion-cyclotron resonance has already taken place 
at lower altitudes. Such a preferential heating is still at work 
between 57" and 102" (cf. Fig. ITTb. If the net heating experi- 
enced by the lowest q/m ion species, as well as their tempera- 
tures, appear high, one should keep in mind, though, that these 
are minor ions. Therefore, this heating does not require much 
energy, nothing comparable to what is needed to appreciably 
heat the protons. And yet, even if we have no direct measure- 
ment of the proton temperature here, one can speculate that 
it is not very different from those of the ions species having 
the higher q/m ratios. One can firmly suppose from the discus- 
sion in the previous paragraph, indeed, that the potential cy- 
clotron heating for the high q/m species is low. Then, it is not 
likely that large temperature differences with the protons can 
be maintained against collisions. If the proton temperature is 
as high as those of the colder ion species observed, this is far 
more than the electro n temperature usual l y observed in coronal 



holes 8 x 10 3 K, Habbal et al. 1993; Wil helmet all 1 19981 



David et al.lll998l) . Then, one has to consider that the protons 
may have been heated in comparison with the electrons, again 
at altitudes below 57". Investigating the helium or hydrogen 
temperatures and their variations with the altitude would be 
a big help for getting more information in the q/m analysis, 
but the associated lines in the SUMER spectrum suffer from a 
strong stray-light blending. 
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Fig. 11. Difference in temperature between 57" and 102", once the non-thermal contribution is removed (left panel). The right 
panel corresponds to the case where we suppose that there is no variation in non-thermal velocity, only a temperature one. This 
case is equivalent to the case ^ = 0, i.e. when no non-thermal velocity is considered at all. The fitting functions are discussed in 
Sect. 15.21 The slopes of the linear fits are -1.9 + 3.1 and -2.4 ± 2.1, respectively, in the left and the right panels (in 10 6 K per 
unit of normalized q/tri). 
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Fig. 12. Temperatures as functions of the charge-to-mass ratio q/m, 57" (left panel) and 102" (right panel) above a polar coronal 
hole (data set 1; the non-thermal velocity has been removed from the line widths). The interest of making a linear fit is discussed 
in Sect. 1331 



5.3. Derived non-thermal velocity and its effect on the 
signature 

The signatures of preferential heating are sensitive to the de- 
rived non-thermal velocity. In particular, any variation in width 
is distributed among variations in temperature and in non- 
thermal velocity: 



A 2 k 

Scr=—-(-6T + t60- 
L(tc l m 



(7) 



Even if there is little variation in £ between 57" and 102" 
(2 kms~'), the actual value of £ has a strong effect on the in- 
terpretation of the line-width variation. The value of %(r) was 
derived under flux conservation condition, with no damping oc- 
curring for the Alfven waves. We have already said that it was 
not possible to verify this assumption, because the Mg x 625 A 
line width keeps increasing, so a possible decrease in £ may be 
hidden by an increase in temperature. If 5%(r) starts to decrease, 
in reality, above some position r or even becomes negative, i.e. 



the damping is greater than the amplitude increase due to the 
flux conservation, then the line-width growth only has to be at- 
tributed to a growth of temperature. But this situation may only 
amplify the preferential heating trend, both in absolute temper- 
atures and in differences in temperature, as shown in the right 
panel of Fig.[TT] Note that the case 5% — is equivalent to the 
case £(ro) = kms -1 (cf. Eq|7]). In contrast, if £(r) is larger 
than what we find, considering for example that Alfven waves 
are not responsible for the apparent non-thermal velocity, hence 
that ^ is not constrained by the flux conservation, then the pref- 
erential heating trends are weakened. 

The non-thermal velocities that we got at low altitude in 
the coronal hole are within the lower limit of previously pub- 
lished results. However, most of those results were obtained by 
using the formation temperature hypothesis, which underesti- 
mates the ion temperatures. 

A value of 25 or 30 kms -1 for the alfvenic fluctuations is 
typically used as a boundary condition at the coronal base in 
numerical simulations of heating and acceleration of the fast 
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Table 3. Line widths measured in data set 1 (coronal hole) at two positions above the solar limb, and the associated ion tempera- 
tures. 









57" 


(ne = 4.6 


x 10 7 cnT J ,£ = 


15 kms~' 


) 




102" 


n e = 2.7 x 10 7 cm' 




17 kms~ 


) 


Line 3 


i h 

q/m b 


Observed width c 




Width d 


T 




Observed width 




Width d 


V 










(A) 




(A) 


(10 6 


K) 




(A) 




(A) 


(10 6 


K) 


Fex 1028.02 


0.16 





1111 


± 0.004 





1102 ±0.004 


5.96 ± 


0.70 





1161 ± 


0.004 





1184 ± 


0.004 


6.76 ± 


0.80 


Fex 1463.50 


0.16 





1434 


± 0.004 





1435 ±0.013 


4.87 ± 


1.15 


















Fexi 1467.08 


0.18 





1310 


± 0.004 





1272 ± 0.007 


3.64 ± 


0.62 





1376 ± 


0.004 





1318 ± 


0.011 


3.73 ± 


1.00 


Ar vm 700.26* 


0.18 





1158 


± 0.004 





1351 ±0.002 


3.35 ± 


0.36 





1222 ± 


0.005 





1411 ± 


0.004 


3.53 ± 


0.58 


Ar vm 713.82* 


0.18 





1291 


± 0.004 






2.88 ± 


0.34 


















Fexn 1242.00 


0.20 





1035 


± 0.004 





1039 ± 0.004 


3.34 ± 


0.50 





1132 ± 


0.004 





1128 ± 


0.004 


3.83 ± 


0.59 


Fexn 1349.36 


0.20 





1071 


± 0.004 





1078 ± 0.004 


2.98 ± 


0.45 





1176 ± 


0.004 





1177 ± 


0.004 


3.46 ± 


0.55 


Cax 574.00* 


0.22 





1111 


± 0.004 





1118 ±0.002 


3.43 ± 


0.42 





1137 ± 


0.004 





1185 ± 


0.002 


3.76 ± 


0.49 


Sivm 1445.75 


0.25 





1605 


± 0.004 






3.24 ± 


0.27 





1674 ± 


0.004 








3.44 ± 


0.31 


Sx 1196.20 


0.28 





1250 


± 0.004 






3.23 ± 


0.35 





1266 ± 


0.004 








3.19 ± 


0.39 


Nv 1238.81 


0.29 





1532 


± 0.004 





1709 ± 0.004 


2.60 ± 


0.18 





1448 ± 


0.004 





1804 ± 


0.004 


2.86 ± 


0.21 


Ovi 1031.93 


0.31 





1426 


± 0.004 






2.98 ± 


0.24 





1406 ± 


0.004 








2.82 ± 


0.26 


Mgrx 706.05* 


0.33 





1669 


± 0.004 





1669 ± 0.002 


3.22 ± 


0.25 





1775 ± 


0.004 





1795 ± 


0.002 


3.67 ± 


0.29 


Na ix 681.72* 


0.35 





1561 


± 0.004 





1587 ± 0.002 


2.94 ± 


0.24 





1535 ± 


0.004 





1637 ± 


0.002 


3.06 ± 


0.27 


Si xi 580.91* 


0.36 





1269 


± 0.004 





1271 ± 0.002 


3.13 ± 


0.32 





1251 ± 


0.004 





1204 ± 


0.002 


2.66 ± 


0.34 


Sixi604.15* f 


0.36 





1229 


± 0.005 






2.66 ± 


0.31 





1313 ± 


0.005 








2.97 ± 


0.40 


Mgx 624.94* 


0.37 





1477 


± 0.004 





1485 ± 0.002 


3.26 ± 


0.28 





1513 ± 


0.004 





1524 ± 


0.002 


3.35 ± 


0.31 



a Wavelengths in angstroms are taken from lFeldman et all dl997h . Lines marked with an asterisk are observed in second order. 
b Ion charge-to-mass ratio. 

c Half-width at 1/ ^[e (Gaussian). Instrumental width is removed, but spectra are not corrected from the instrumental stray light. For lines 

seen in second order, the widths are twice as large as what would be measured in first order (see Sect. l5.2t . 
d Half-width at 1/ ^fe when spectra are corrected from the stray light, with the instrumental width removed. 

e Ion temperature deduced from the corrected line width, with the non-thermal velocity removed. When no correction was necessary or 
possible, the observed width was used. 

This line is blended with an Fe vn line, as given bv lEkberg & FeldmarJ l2003h . See comments in Sect. 13. II 



solar wind (see references in Sect. Q]). However, this value is 
only encountered above 200" in our results. There is then a 
lack of energy at the coronal base to reproduce the speed and 
temperature of the fast solar wind by using these models. 



6. Conclusions 

Alfven waves and ion-cyclotron resonance, i.e. absorption of 
the high frequencies, are commonly used in theoretical mod- 
els to provide the corona with heating and/or the fast solar 
wind with acceleration. Signatures of ion-cyclotron heating 
have been reported in the upper corona (> 2.5 R Q ) and in the so- 
lar wind. It appears important, then, to investigate ion-cyclotron 
resonance closer to the Sun, where the coronal line widths of- 
fer the opportunity to study both the temperature of different 
ion species and the amplitude of possible Alfven waves. This 
is of course a challenging task as these quantities are merged 
into one observable, the line width. Most of the time, methods 
were used that excluded a priori the possibility of analysing 
both contributions simultaneously. Our analysis, on the con- 
trary, achieves a synthesis of different approaches that can be 
found separately in the literature, to both test and constrain the 
presence of Alfven waves and the preferential heating. We con- 
strained the former by the gradient of the Mgx 625 A line 
width, while the latter appears from the gradient of ion tem- 



peratures and from the distinct ion temperatures measured at a 
given altitude. 

Our analysis rests on the following hypotheses: i) the coro- 
nal lines broadening is solely due to the thermal and non- 
thermal Doppler effects through Gaussian contributions, ii) The 
non-thermal velocity ^ comes from the integration of fluid mo- 
tions of the coronal plasma on both the observed solid angle 
and the exposure time, and is identical for all ion species at 
a given altitude, iii) These motions are due to the presence of 
Alfven waves; if there is no damping, their amplitude increases 
with the altitude to conserve the energy flux in a density- 
stratified corona, iv) The increase in width of the Mg x 625 A 
line (the highest q/m), at low altitude, results only from un- 
damped Alfven waves, i.e. Ar = 0, so that a boundary con- 
straint can be set for £ using a density diagnostic (Si vm line 
ratio) and the energy flux conservation, v) The Alfven waves 
stay undamped on the entire observed spatial range, so that %(r) 
can be derived. 

Several results must be emphasised. The decrease in the 
coronal line widths, occurring around 1.2 R for Mgx for ex- 
ample, is essentially an effect of the instrumental stray light 
contamination that disappear when we remove the predicted 
stray-light spectrum. We cannot rule out, though, the existence 
of a plateau at higher altitudes than we observed. This result is 
in accordance with some theoretical ones tha t find it diffic ult to 
damp the Alfven waves in the lower corona dParkerlll99ll) . 
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The derived non-thermal velocity is about 15 km s _I at 60", 
while a value of 25 kms -1 is only reached at 200" above the 
limb. This means a widely lower wave energy than that com- 
monly used in numerical simulations of heating and accelera- 
tion of the fast solar wind. Consequently, one has to find an- 
other form of energy, or of dissipation process, to replace or 
at least complement the processes usually presented in models 
based on alfvenic fluctuations. 

Two kinds of signatures of preferential heating are found: 
ion species having the lowest q/m ratios are the hottest ones 
at a given altitude and experience more heating as the altitude 
increases from 57" to 102". If some damping of the Alfven 
waves occurred, i.e. if hypothesis v) was wrong, then these 
trends would be stronger. Yet, this is neither a proof by itself of 
cyclotron resonance occurring in the corona, nor does it imply 
that this is a dominant energy source for acceleration of the fast 
solar wind or the heating of the corona: the energetics have to 
be analysed in more depth. Minor ions can only help to trace the 
cyclotron resonance that may also concern the protons. The ion 
temperatures and probably the proton one are higher than the 
electron temperature usually measured in coronal holes, which 
imply some preferential heating that must have taken place be- 
low the lowest altitude where we observed, i.e. 57". 

Some facts may possibly weaken our analysis. First, the 
width of some lines starts to decrease at the lowest altitude we 
observed. This contradicts our interpretation in terms of a com- 
mon increase in the wave amplitude, but this can be due to 
another contribution in the line width, influencing these lines 
more. Second, the lines associated with the lowest q/m ratios 
(Fe x and Fe xi, particularly) are highly contaminated with stray 
light, so that it is difficult to follow their behaviour at higher al- 
titudes. We are also aware of the large error bars in our derived 
quantities. Therefore, our results have to be confirmed by addi- 
tional observations and further analyses. 

This work also emphasises some needs for future instru- 
mentation, especially in terms of spectrometers. The study of 
line widths emitted in the low corona requires a small instru- 
mental width as compared to the intrinsic width. It is important 
to develop a better rejection of the light coming from the so- 
lar disc, for off-limb observations. Furthermore, a solar probe 
could achieve important in situ observations to complement the 
remote-sensing observations, and even calibrate their interpre- 
tation. 
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